Dynamic Modes of Microcapsules in Steady Shear Flow: Effects of Bending and Shear 

Elasticities 
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The dynamics of microcapsules in steady shear flow was studied using a theoretical approach 
based on three variables: The Taylor deformation parameter qd, the inclination angle 9, and the 
phase angle (j) of the membrane rotation. It is found that the dynamic phase diagram shows a 
remarkable change with an increase in the ratio of the membrane shear and bending elasticities. 
A fluid vesicle (no shear elasticity) exhibits three dynamic modes: (i) Tank-treading (TT) at low 
viscosity lyin of internal fluid (q:d and 6 relaxes to constant values), (ii) Tumbling (TB) at high r^n 
{6 rotates), and (iii) Swinging (SW) at middle r/in and high shear rate 7 {9 oscillates). All of three 
modes are accompanied by a membrane {(f>) rotation. For microcapsules with low shear elasticity, 
the TB phase with no (j) rotation and the coexistence phase of SW and TB motions are induced by 
the energy barrier of </> rotation. Synchronization of (p rotation with TB rotation or SW oscillation 
occurs with integer ratios of rotational frequencies. At high shear elasticity, where a saddle point in 
the energy potential disappears, intermediate phases vanish, and either 4> or 6 rotation occurs. This 
phase behavior agrees with recent simulation results of microcapsules with low bending elasticity. 

PACS numbers: 87.16.D-, 05.45.-a, 82.40.Bj 



I. INTRODUCTION 

Soft deformable objects such as liquid droplets [H, [2], 
lipid vesicles i-iiLred blood cells (RBCs) [Il-lp, and 
synthetic capsules p3 - [5l| exhibit rich behaviors in flows. 
In recent years, they have received growing attention ex- 
perimentally, theoretically, and numerically. These ob- 
jects are filled with liquids and have a constant internal 
volume V . The lipid vesicles, RBCs, and synthetic cap- 
sules are wrapped by lipid or polymer membranes and 
resist breaking up under much larger shear than liquid 
droplets. The variation of material properties of these 
objects can lead to different dynamic modes with mor- 
phological changes. 

Since microcapsules can carry substances and control 
their release, they are widely used in various fields in- 
cluding food, printing inks, pharmaceutics (drug carrier 
for drug delivery system), agrochemicals, and cosmet- 
ics. Various kinds of microcapsules are synthesized with 
different sizes and membrane thicknesses h 42, 43, 52.^ 
Is^ l- The deformability of microcapsules can be varied 
widely by controlling the membrane elastic properties 
and capsule shapes. When the membrane is made of 
linear-elastic material, two-dimensional Young's modulus 
l^nb = 4:fiKp^/{fi + KfC) = Yh and the bending modulus 
K = Yh^ /\2{\—v'^), where /i and Kx are two-dimensional 
shear and area compression moduli, respectively, with 
three-dimensional Young's modulus Y and Poisson ratio 
V [ssl ]. Thus, thicker membranes become more difficult 
to bend than to stretch {k/Y^y> oc h?)- 

RBCs have a composite membrane made of lipid mem- 
brane and attached spectrin network. The lipid mem- 
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brane is in a fluid phase and is almost incompressible. 
The spectrin network adds shear elasticity to the mem- 
brane. The bending and shear elasticities of RBCs are 
estimated by micropipette aspiration p6l [56^ and opti- 
cal tweezers experiments [57i - i59l] : k ~ 2 x 10~^^J and 
/i ~ 6 X lO^^N/m. 

The relative importance of bending and shear elas- 
ticities is determined by the dimensionless Foppl-von 
Karman number 7 — Y^\,R^Jk, where i?o is the mean 
radius of capsules or cells [60|- The buckling transi- 
tion of capsules from spherical to icosahedral shape oc- 
curs at 7b ~ 150 60]. Capsules with a linear elas- 
tic membrane have large values of 7 > 100 because 
7 = 12(1 — i^^)(i?o/^)^- Capsules with a composite mem- 
brane can also have small 7. In this paper, we investigate 
microcapsules with wide ranges of 7. 

Recently, the dynamics of lipid vesicles in steady shear 
flow was intensively investigated (3l-[2^. A lipid vesicle 
can be considered as a microcapsule in the small limit 
of the shear modulus /i — 0. An isolated fluid vesi- 
cle exhibits three types of dynamic modes in a steady 
shear flow, with flow velocity v = jyex, where is 
a unit vector in the flow direction. When the viscos- 
ity of internal fluid yyin and membrane viscosity Tymb are 
low, the membrane exhibit a tank-treading (TT) rota- 
tion with a stationary shape and a constant inclination 
angle 6 > [see Fig. [Hb)]. At high or rymb, the vesi- 
cle exhibit a tumbling (TB) motion, where 9 rotates [see 
Fig. [IJc)]. Around the TT-TB transition viscosity with 
hi gh shear rate 7, a swinging (SW) motion appears 0- 
[qI. Tisl. [Tsl - fTTi [2l| ■ where 6 and the vesicle shape oscillate 
see Fig. Hfd)]- This motion is also called trembling 0- 
m, Il3 or vacillating-breathing [H [11]. These three 
types of motion can be understood by the perturbation 
theories for quasi-spherical vesicles [15l - [l7| or a general- 
ized Keller-Skalak (KS) theory for deformable ellipsoidal 



vesicles [2T|. 

Experimentally, synthetic capsules and RBCs show the 
oscillation (swinging) of their lengths and 9 [l^, H^, lisj 
during TT motion and RBCs also transit from TB to 
TT with increasing 7 (28l . |6]| . Recently, this dynam- 
ics was explained by the KS theory with the addition 
of an energy barrier for the TT rotation caused by the 
membrane shear elasticity [1^, H^. More recently, we 
phenomenologically extended this theory [sl] to include 
the shape deformation of RBCs [41']. Fix-shaped and de- 
formable RBCs give very similar phase behavior. At low 
7, RBCs show TB [0) rotation without TT (<^) membrane 
rotation. In the middle ranges of shear rate 7, synchro- 
nized phases of the 9 and rotations with integer ratios 
of the rotation frequencies as well as intermittent rota- 
tions |36|] appear [4l|. The TT-TB transition was also 
obtained by recent simulations of microcapsules jSjiGj . 
However, the above intermediate phase (synchronized or 
intermittent rotations) was not observed in their simula- 
tions. Kessler et al. contested that intermittent rotation 
is an artifact of the theoretical model [13, E^l . As we de- 
scribe below, the disappearance of the intermediate phase 
is found to be caused by low bending elasticity in their 
simulations. In this paper, we systematically investigate 
the dynamics of microcapsules in steady shear flow, and 
show that their phase behaviors are qualitatively changed 
from the bending-elasticity dominant condition to the 
shear-elasticity dominant condition. 

Note that some confusion is caused by the different def- 
initions of a tumbling state. For fluid vesicles, tumbling 
typically means the rotation of inclination angle accom- 
panied by the rotation of the phase angle 0. On the other 
hand, for RBCs and capsules, it typically means the 9 
rotation accompanied by the oscillation (not rotation). 
Thus, in the definition for fluid vesicles, synchronized or 
intermittent rotations of 9 and are also called "tum- 
bling". In this paper, we call them "TB with no TT" 
and "TB with TT" to distinguish their modes. 

The dynamics of microcapsules is described by several 
dimensionless quantities. Here, dimensionless quantities 
are denoted by a superscript *. The relative ratio of the 
volume V and surface area S is characterized by the re- 
duced volume V* = 3F/(47rS')3/2 = {Ry/Rsf = (1 + 
As/4^)-3/2, where i?v = {W/At:^/^, Rs = (5/47r)i/2, 
and the excess area As = S/Ry — 47r. The relative vis- 
cosity of the inside fluid and membrane are rj*^ = ?7in/'?o 
and r/*^^ — rjaih/voRs, where rjo is the viscosity of the 
outside fluid. The bending rigidity k is used as an en- 
ergy unit, and the shape relaxation time of the fluid vesi- 
cles with bending rigidity k is given by r = tjoRc^/k (for 
77i*j — 1). This time is used to define a reduced shear rate 
7* — 7T. Since the Reynolds number is low, Re< 1 in 
typical experimental conditions, the effects of the inertia 
are neglected. 

The theoretical models of fluid vesicles and microcap- 
sules are explained in Sec. [Hi The results are described 
in Sec. IIIII First, the dynamics of fluid vesicles — 0) 
is briefly explained in Sec. IIII Al and then dynamics 
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FIG. 1: Schematic representation of (a) a capsule and its 
dynamic modes: (b) Tank-treading (TT), (c) Tumbling (TB), 
and (d) Swinging (SW) modes, (a) The inclination angle d 
and phase angle (/> are depicted. The closed and open circles 
are the original and current positions of a membrane segment, 
respectively, (b) For an elastic capsule, a capsule shape and d 
oscillate, while they are kept at fixed values for a fluid vesicle. 



with low and high shear elasticity [i is described in Sees. 
IIII BfjHlDl The mechanism of the dynamic modes and 
experimental conditions are discussed in Sec. IIVI The 
summary is given in Sec. |Vl 



II. GENERALIZED KELLER-SKALAK THEORY 
A. Fluid vesicle 

Keller and Skalak [sO] analytically derived the equation 
of motion of the inclination angle 9 for vesicles or elastic 
capsules based on Jeffery's theory |62|. The KS theory 
assumes (i) a fixed ellipsoidal shape of vesicles and (ii) a 
simple velocity field on the membrane. It gives very good 
predictions for vesicle dynamics at low shear rates. The 
KS theory is phenomenologically extended to include the 
shape deformation of fluid vesicles on the basis of the per- 
turbation theory [l2,[l3,[lB,[i7| of quasi-spherical vesicles 
in Ref . [2l| . The shape deformation is represented by the 
Taylor deformation parameter an = {Li — L2)/{Li+L2), 
where ^1,^2 are the principal lengths on the vorticity 
{xy) plane, and L3 is in the vorticity (z) direction. The 
differential equations of q;d and 9 are given by 



jdt 



= <l- 



-f = i{-l + /oAcos(2.)} 
d0 _ cos(2e') 



jdt 




2/i{l + /2(??r„-l) + /2M;b} 



where ^0 = 45/87r(32 + 23?7*^ 16r4i,)F* and Ai = 
30/(32 + 237/i*„ + 16<b) III- The first and second 
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fixed ellipsoidal shapes [S^. Equation ^ is modified as 



FIG. 2: Free-energy potential function of the Taylor defor- 
mation parameter qd at </!> = and (p — n/2 [F{—ao,0) ~ 
F(aD,7r/2)]. Solid and long-dashed lines represent fee = 
and kc — 10 with fes — 0, respectively. Short-dashed line rep- 
resents kc — 10 with fcs = 1. Dashed-doted line represents the 
potential of the shear elasticity estimated by the perturbation 
of quasi-spherical capsules with fiRg/K = 1. 



terms in the last parentheses of Eq. ((TJ represent the 
forces of the membrane bending elasticity and the exter- 
nal shear stress, respectively. The free energy is given by 
F = Fcv = (k/2) / (Ci + C2fdS, where Ci and C2 are 
the principal curvatures at each point of the membrane. 

In Eq. ([5]), the factors /o, /i, /2, and /s are function 
of the vesicle shape (ao, V*) (see Appendix A). Because 
the factors and F are independent of the phase angle 4> 
for fluid vesicles, Eq. ([2]) can be simply written as 



dO 



1 



{-l + Bcos(26')}. 



(4) 



The factor _B is a function of od, V* ^ rj^^, and fy^i,. 

At rj*^ = 1, the lipid vesicle shows only the TT motion 
at any V* and 7*. The TT steady angle 9o predicted 
by KS theory [Eq. ^] gives very goo d agreement with 
simulations and experiments 0, 

iinillfllll, US- On 

the other hand, the perturbation theories give TB mo- 
tion even at 7^;*^^ = 1 for 0.8, while they give good 
predictions for V* > 0.9 [2j]. Thus, the perturbation 
theories should not be applied to the dynamics of vesi- 
cles and capsules for V* < 0.8. 



B. Microcapsule 

Synthetic capsules and RBCs have membrane shear 
elasticity /i, which modifies their free-energy landscape: 
F — Fcv + F^. The energy of the shear elasticity 
depends on the capsule shape and phase angle 0, while 
Fcv only depends on the capsule shape. The angle </) = 
is the original position in the absence of flow. During 
(j) (TT) rotation, the membrane is locally elongated and 
it generates a free-energy barrier at (f) = 7r/2. Recently, 
Skotheim and Secomb extended the KS theory to take 
into account this free-energy barrier for capsules with 



d(t> _ {co/j*)dF*/d(j} + cos{20) 

^ " ' 2/i{i + MvL - 1) + /2/3Cb} 



(5) 



where cq = S/a/Svr/iF* [36]. 

In this paper, we employ Eqs. ([1]), (HI, and ([S]) to in- 
vestigate the dynamics of the elastic capsule with pro- 
late shape at V* = 0.9 and r]*^^^ — 0. These equa- 
tions are numerically integrated using the fourth-order 
Runge-Kutta method. The curvature energy i^cv(aD) is 
numerically calculated for the ellipsoidal vesicle. The 
prolate (Li > L2 = L3) and oblate (Li = L2 > L3 and 
ttD = 0) shapes are energy minimum and maximum, re- 
spectively (see Fig. (2). Here, a fit function is employed to 
avoid the artifacts of nonsmooth functions: dF*^/ daj) = 



30aD 



820a^ 



ckD + 1.7/yaf^ with 



l.T/^ag^-^ - ckD + i-'/V": 
Q,max _ Q 35g Here, the relative energy to the free- 
energy minimum (prolate shape) is considered. The force 
—dFcv/dau diverges in the limit of maximum exten- 
sion, ttD — >■ a^^^{V*). For shear elasticity, the leading- 
order energy function is employed in Sees. IIII Bl and 
ImCl Ff,{ao,(t>)* = (fcc/2){(aD - aofcos''{^) + {au + 
ao)^ sin^ ((/))}, where ao = 0.325 is an of the prolate 
shape. Because compression and 7r/2 rotation can gen- 
erate the same capsule state, the free energy should keep 
the relation _F(— aD,0) = F(aD, + 7r/2). The shear 
elastic energy is estimated by the perturbation of quasi- 
spherical capsules; see Appendix B. The harmonic po- 
tential i^*(Q;D,0) with fee — 4/ii?g/K gives good agree- 
ment with the perturbation result; see Fig. [2l In Sec. 
IIIID} in order to check dependence on the potential 
shape, a higher-order term ^3 is added to the potential: 



F3* = (A:efc3/3){(aD-«o)' 



^(0)-(aD + ao)^ 



III. RESULTS 



A. Fluid vesicle 



First, we briefly describe the dynamics of lipid vesi- 
cles (fee — 0) in steady shear flow. Figure ^a) shows 
its phase diagram. At low viscosity ratio r]*^, the vesicle 
shows a TT motion, where an and remain constant. 
The factor B in Eq.^ is more than unity, and a sta- 
ble fixed point 9 = do = 0.5 arccos(l/B) and unstable 
fixed point = —Oq exist. As rj*^ increases, B and 6*0 
decrease. At the critical viscosity = = 5.27, B 
becomes unity and the two fixed points 6 — ±0o merge 
at 6* = (saddle- node bifurcation). At 77*^^ > 77* with 
low shear rate 7* < 1, no fixed point exists so 9 rotates 
(TB with TT); see Figs. [SJb-d). This TT-TB transition 
is quantitatively predicted by the original KS theory for 
fix-shaped vesicles. 

On the other hand, at higher shear 7* > 1, the shape 
deformation of the vesicles is not negligible. A SW phase 
appears between the TT and TB phases. In the SW 
phase, 9 and a-Q oscillate; see dashed lines in Figs. EJb- 
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FIG. 3: (Color online) Dynamics of fluid vesicles (fcc = 0). 
(a) Phase diagram (??*„, 7*). Solid (red and brown) lines rep- 
resent phase boundaries between TB, TT, and SW. For the 
right region of the dashed line, the vesicle passes through 
or remains at ^ < 0. (b-d) Time development of od and 
Q. Solid and dashed lines represent the TB and SW mo- 
tions at (?7in,7*) = (6,1) and (6,1.2), respectively. Dashed- 
dotted lines and closed circle represent the TT motion at 

(^in,7*) = (5,1). 



d). These oscillations are generated by the shape defor- 
mation as follows 21]: In Eq. ((T|), the shear force de- 
pends on Q as sin(26'), so the shear increases an (elonga- 
tion) for < 6* < 7r/2, but decreases ao (shrinkage) for 
— 7r/2 < < 0. A prolate vesicle starts rotation with 
_B < 1 as in the TB phase and then it shrinks to a more 
spherical shape at 6* < 0, which has greater B. When 
B becomes greater than 1/ cos(26'), the right-hand side 
of Eq.([l]) changes its sign and the angle d increases. At 
> 0, the vesicle elongates back to the prolate shape. 
The TB-SW transition occurs when the trajectory 
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FIG. 4: (Color online) Phase diagram (??*„, 7*) of microcap- 
sules at fee = 1. (a) Enlarged figure for the region of synchro- 
nized rotation, (b) Region at high 77*^. (c) Region at low 7*. 
Solid lines represent the phase boundaries. In the right region 
of dashed line, the vesicle passes through or remains at S < 0. 
Numbers represent the rotation frequencies /^ot : ^to^■ 



crosses an — 0. Between 77*^ = 77* (dashed line in Fig. 
[3l[a)) and the SW phase, the TT phase has negative 
steady angle 6*0 < 0. For fluid vesicles, the phase angle 
always shows rotation except for the limit 77*^^ — >■ 00. 
At ?7j*jj — > 00, Eq. becomes Jeffery's equation fG^] for 
solid ellipsoidal objects, where is kept in the original 
position. These results show good agreement with recent 
experiments of lipid vesicles 



B. Microcapsule with low shear elasticity 

When the membrane shear elasticity is added, the dy- 
namics is changed from the above dynamics of fluid vesi- 
cles. First, we consider the dynamics of microcapsules 
at fco = 1, where shear elasticity is relatively lower than 
the bending elasticity. The energy barrier for rotation 
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FIG. 5: (Color online) Dynamic mode transition from TB to 
TT phase at fee = 1 and rji^ = 5. (a) Rotation frequencies 
/rot of 9 (solid line) and <f) (dashed line), (b) Trajectories of 
TB with fL ■■ fL = 1 : and TT with f^,, : ft, =0:1 
at 7* = 0.025 (dashed line) and 0.07 (solid line), respectively, 
(c) Trajectory of TB accompanied by TT rotation with 
/rot : fL = 1 : 2 at 7* = 0.05. 



is F* — 0.21 ai (f) — 7r/2, which is five times lower than 
the energy barrier F* = 0.97 for the capsule shrinkage 
at ttD = 0. Because of this energy scale difference, the 
phase diagram can be separated into two regions: high 
shear region 7* > 0.8 shown in Figs. IDJa) and (b), and 
low shear region 7* < 0.1 shown in Fig. SJc). 
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FIG. 6: (Color online) Dynamic mode transition from TB to 
SW phase at fce = 1 and rj*-^ = 6. (a) Rotation frequencies 
/rot of 9 (solid line), (j) (dashed line), and swinging (dashed- 
dotted line), (b) Trajectory at 7* = 0.85: TB accompanied 
by TT (j) rotation, (c) Trajectory at 7* — 0.95: Coexistence 
of TB and SW motions, (d) Trajectory at 7* = 1.05: SW 
motion, (e) Return map at 7* = 0.95 for the angle <j> when 
trajectory crosses 9 — at qd > 0.2. Dashed line represents 

0n + l = (pn- 



1. TB-TT transition 

At low shear rate 7* < 0.1, the shape deformation 
is very small so microcapsules keep their prolate shape, 
i.e. ckD — OLQ. Thus, the generalized KS theory for fix- 
shaped capsules proposed by Skotheim and Secomb [SG^ 
(Eq. ([2|) and Eq. ([5]) with constant an) predicts the 
phase behavior in this region very well. At very low shear 
rate 7* < 0.02, the capsules show TB motion with no 
TT unlike fluid vesicles; see dashed line in Fig. [SJb). 
The energy barrier prevents rotation. Without TT 
rotation, the whole capsule has to rotate (tumble) due to 
the shear torque applied on the capsule surface. 

As 7* increases, (p starts to rotate and then the TT 
frequency fL gradually increases and the TB frequency 
fL decreases; see Fig. [SJa). Here, an angle change of 
TT is counted as one rotation. In this coexistence range 
of two rotations (TB with TT) , the synchronization of 
and 9 can occur with the integer ratios of fL and fL- 
This infinite number of synchronization phases is called 
Devil's staircase [g^l- The trajectory of the synchronized 



rotations with fL ■ fL = 1:2 are shown in Fig. [Sfc). 

At sufliciently large 7* > 7*^ (7*^ = 0.068 at r]^^ = 5), 9 
stops rotation and the TT phase appears. With the fre- 
quency of (f) rotation, 9 and ao oscillate (swing) around 
the steady values for fluid vesicles (fco = 0). Asjy* in- 
creases, the oscillation amplitudes decrease 0, |4l|. In 
this paper, we call this phase TT following the phase di- 
agram of fluid vesicles (it is also called SW [2^, 13, HI] 
or OS in the previous papers). As 7* — >■ 7^*^, the 
mean angle (9) approaches to null and the amplitude 
increases and then passes through < during this 
oscillation; see solid line for 7* = 0.07 in Fig. El^b). 
This crossing region appears between the solid (red) and 
dashed lines in Fig. lUJc). This dashed line does not rep- 
resent a phase boundary. In crossing this dashed line, 
the averages and amplitudes of 9 and ao show smooth 
change 14l|. No discrete changes of trajectories are seen. 
This trajectory crossing 9 < resembles the one in the 
SW phase of fluid vesicles. However, they are caused by 
a different mechanism. We call them type 1 and type 11 
(swinging) oscillations. Type I: Oscillation induced by 
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FIG. 8: (Color online) Trajectories of microcapsules at fco = 1 
and 7* = 1.2. (a) Dashed and solid lines represent trajectories 
in the TT phases (type II oscillation) at 77;*^ = 5.748 and in 
the SW phase (unsynchronized type I and II oscillations) at 
r/in — 5.7486, respectively, (b) Unsynchronized oscillations 
with larger amplitude of the type I oscillation in the SW phase 
a-t J^'n = 5.75. (c-d) Synchronized rotations in the SW phase 
at ryfj-j = 6.12 (dashed line) and in the TB phase at 'ql^ = 6.14 
(solid line). 



FIG. 7; (Color online) Dynamic mode transitions from TT, 
SW to TB phase via synchronized region at fee = 1 and 7* = 
1.2: TT for r?*„ < 5.7484, SW for 5.7484 < r)*^ < 6.127, 
SW+TB for 6.154 < 77*, < 6.219, and TB with TT for 6.127 < 
r^fn < 6.154 and 77;'^ > 6.219. Maximum and minimum of (a) 
qd and (b) 9. The region around transition between TT and 
SW is enlarged in the inset in (a), (c) Rotation frequencies 
/rot of 9 (solid line), (j> (dashed line), and swinging t/j (dashed- 
dotted line). 



the shape deformation. For the oscillation of fluid vesi- 
cles, the large shape deformation is essential as explained 
in Sec. IIII Al Type II: Oscillation induced by the energy 
barrier to cf) rotation. This oscillation can be reproduced 
with a fix-shaped approximation including the passage 
through the negative angle 9. The resulting oscillation 
of the (j) rotational velocity causes 9 and an oscillations. 
These two types of oscillations cannot be distinguished 
only by trajectory shapes. 



2. SW phase 

The microcapsule at fee = 1 has a SW phase around 
the region of the SW phase at fee = 0; compare Figs. (SJa) 
anddfb). In the SW phase, the type I oscillation occurs 
in addition to the type II oscillation. It is found that 
this coexistence makes the transition from TB with TT to 
SW a gradual change; see Fig. [B] The SW ■0 frequency is 
calculated as a frequency of crossing 9 = 9m\d upwards at 



Q!D < Q!mid in Fig- HI where (6'mid, amid) is a point inside 
the SW trajectory, and (0mid,amid) = (0,0.2) is chosen 
here. The type II oscillation generates the fluctuations 
of the TB or SW trajectories, when the oscillations and 
(j) rotation are not synchronized; see return map in Fig. 
[S{e). As 7* increases, an unsynchronized (intermittent) 
SW oscillation starts, and then its frequency increases 
while the total of the TB and SW frequencies is almost 
constant. 

The TT-SW transition is slightly modified from fluid 
vesicles. Both TT and SW phases are accompanied by 
the type II oscillation. The transition point is detected 
as a discrete point of the slope of the amplitudes of q;d or 
9 with respect to 77*^^; see the region around the transition 
point 77i*„ = 5.7484 in Figs. Elja) and (b). Two types of 
oscillations without synchronization are shown in Figs. 
Ha) and (b). 

Synchronization of type I and type II oscillations can 
occur with integer ratio of the oscillation frequencies. 
The synchronization with /^^^ : ff^^ : ff^^. =0:2:1 has 
the widest region (denoted by 'Sync SW' in Fig. lU^a)). 
Two (p rotations occur every 9 oscillation in this region; 
see dashed fines in Figs. HJc) and (d). The TB-SW tran- 
sition in synchronized region (77*^^ = 6.127 at 7* = 1.2) 
is a discrete transition like that of the fluid vesicles; The 
limit cycle with f^^,. : ff^^. : /,tt = : 2 : 1 becomes 

/rot • /rot ■ /rot = 1:1:0. At the boundaries of the 
synchronized region (77*^^ = 6.094 and 6.154 at 7* = 1.2), 
the trajectories show discrete changes; see Figs. [Tfa- 
c). Other synchronized regions with different ratios of 
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FIG. 9: (Color online) Phase diagram (fee, 7*) of microcap- 
sules at 77in = 5. Top and bottom lines represent the phase 
boundaries of TT {f^^^ : /,tt =0:1) and TB with no TT 
(/rot ■ /rot = 1-0), respectively. Shadowed regions represent 
synchronized phase with : ff^^ = 1 : 2, 1 : 1, and 2 : 1. 




n 
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FIG. 10: (Color online) Phase diagram (J7in,7*) of microcap- 
sules at kc — 5. Lines represent the phase boundaries and the 
region with 61 < as depicted in Fig. U 



/rot ■ /rot seen as kinks or slits of the boundary line 
between 'SW+TB' and 'TB with TT' phases in Fig. |4i;b). 



C. Microcapsule with high shear elasticity 

With increasing fee , the phase diagram is found to show 
remarkable changes. Phase diagrams at rj*^^ = 5 and at 
ke = 5, 10, and 100 are shown in Figs. I9HIT1 At < 5, 
the TT-TB transition shear rates 7* are approximately 
linear to kg, as the height of the energy barrier for 
rotation linearly increases; see Fig. [HI A free-energy 




FIG. 11: (Color online) Phase diagram (77i*n,7*) of microcap- 
sules at (a) fee — 10 and (b) fee = 100. Lines represent the 
phase boundaries and the region with 6 < as depicted in 
Fig.H 

minimum at = 7r/2 (or an < at = 0) disappears 
at kc > k^ — 8.2; see Fig. [2l Around this threshold 
kc ~ k^, the phase behavior dramatically changes. At 
kc — 5 (slightly below k^), the TB-TT transition line be- 
comes close to the SW phase boundary as shown in Fig. 
[TUl The larger amplitudes of type II oscillation make the 
coexistence region 'SW -|~TB' wider. The synchronized 
region with f^^^ : ff^^ = 1 : 1 at the TB-SW boundary 
also grows from A77*„ = 0.04 (fco = 1) to A?7*„ = 0.4 
{kc = 5). A slit at i]^^ = 6.9 in Fig. [TUl shows a syn- 
chronized region with f^^^. : ff^^ = 3:2. Narrower slits 
with different ratios of f^^^. : /^^(. are omitted in the phase 
diagram in Fig. 1101 

With a further increase in kc, the region width of 'TB 
with TT' decreases, and then vanishes; see Figs. [S] and 
[TT] At 77*^ = 5, first the regions of unsynchronized (in- 
termittent) rotation are reduced, and the synchronized 
region with f|!^^. : /^^^ = 1:1 becomes dominant. At 
kc ~ 8, this synchronized region only remains for 'TB 
with TT'. At fco > 9, direct transition from 'TB with no 
TT' to TT phase occurs. 

At r;*„ = 100, only two phases, TT with f^^, : ft, = 
: 1 and TB with J^^^ : ff^^. = 1:0, remain; see Fig. 
[TTl b). This phase diagram agrees well with recent sim- 
ulations for microcapsules with high shear elasticity at 
7 = 100 dil and 7 00 [11]. The rotations of 4> and 
9 do not occur simultaneously at any point in the phase 
diagram. During tumbling, the shrinking force caused by 
shear elasticity prevents a passage through <j) = 7r/2. In 
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FIG. 12: (Color online) Viscosity ratio 77;*^ dependence of 
maximum and minimum of (a) qd, (b) 0, and (c) </!> at 
fce = 100. Dashed and solid lines represent 7* = 2 and 4, 
respectively. 



the TT phase, the shear force cx 7sin(26') allows the cj) 
rotation through = tt/2. The transition between these 
two phases occurs when the trajectory reaches au = 
like the SW-TB transition of fluid vesicles; see Fig. [T^l 
In the TT phase, the amplitudes of 9 and au oscillations 
smoothly increase with increasing rj*^; 9 and oscillate 
with the TT frequency {f^^, : ft, : = : 1 : 1) as 
in the TT phase at smaller fcg. At small and large 
type II and type I oscillations may be dominant, respec- 
tively. However, a clear boundary of TT and SW phases 
is not found. The dashed line in Fig. [TlT b'). where the 
minimum of 9 crosses = 0, is not a phase boundary 
either. 



D. Effects of a saddle point of the potential 

As shown above, the phase diagram of microcapsules 
is changed at /cc — ^c, where the saddle point in the free 
energy (the energy minimum at Q!d < for = and 
the energy maximum in the (p rotation for constant q;d) 
vanishes. To check whether the saddle point is essential 
for the phase behavior or it is an accidental coincidence, 
the higher-order term F3 is added to the potential. This 
term restores the energy minimum at ajy < 0, while the 
energy profile is little modified at ao > as shown in 




0.5 kg 1 

FIG. 13: (Color online) Phase diagram (^3,7*) of microcap- 
sules at = 5 and feo = 10. Lines and shadowed regions 
represent the phase boundaries and synchronized phase, re- 
spectively, as depicted in Fig. |9l The coexistence region of 
two limit cycles TT (/4t : ftt =0:1) and TB with TT 
(/rot : /rtt = 1 : 1) is seen around (^3,7*) = (0.3,0.3). 

Fig. [2] Figure [T3l shows that the synchronized and inter- 
mittent regions in the phase diagram reappear with an 
increase in k^. The resulting dynamics are very similar 
to those at low with k^ = 0. Thus, the free-energy 
landscape for cf) rotation is a key feature to determine 
the dynamics of microcapsules. 

IV. DISCUSSION 



Our results revealed that a saddle point of the potential 
changes the phase diagram of microcapsules. For prolate 
capsules with linear-elastic membrane at V* = 0.9, the 
threshold of the phase diagram fee ~ 10 corresponds to 
Foppl-von Karman number 7 = YmhRo/K = 1 10, 
which is lower than those of typical synthetic microcap- 
sules (7 > 100) and RBCs (7 ~ 400). However, the ex- 
istence of a saddle point also depends on capsule shapes 
and buckling instability. A saddle point is observed for 
RBCs by experiments [B^l and simulations [4l[ although 
RBCs have large 7. When an RBC is compressed in the 
X direction (y axis is the symmetric axis of a discocyte) 
with a constraint of (j) — 0, two dimples are buckled and 
new dimples appear on the x axis. These dimples gener- 
ate an energy minimum at </> = (saddle point). There- 
fore, RBCs exhibit synchronized or intermittent rotations 
of0and0[4l|. 

Fischer reported two types of relaxation dynamics for 
RBCs after the shear stops |34t|. After TT motions in 
steady shear flow, the shear is stopped at ~ 7r/2. Then 
the RBCs return to their thermal equilibrium state by 
rotation or shrinkage; see Fig. [T31 During the shrinkage, 
the buckling is observed. For microcapsules, the exis- 
tence of a saddle point of the potential can be checked 
by the observation of this relaxation. If capsules have no 
saddle point, it would relax by continuous shrinkage (not 
by (f) rotation) . Buckling or (p rotation suggests an energy 
barrier for the shrinkage. For synthetic capsules, a sad- 
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^_ rotation _ _ 



^ shrinkage 







FIG. 14: Schematic representation of two relaxation path- 
ways from (aD,0) ~ (ao,7r/2) to (aD,0) = (ao,0) after stop- 
ping shear flow. 



die point can be induced by buckling and nonlinear shear 
elasticity. At large 7 > 150, buckling transition can oc- 
cur for spherical and ellipsoidal capsules [13] • Hysteresis 
at the buckling transition was reported in the simulation 
of microcapsules with 7 — cx) for elongational flow f50| . 
Similar buckling may appear in shear flow. 

We assumed that the major axis of the microcapsule 
is on the vorticity (xy) plane. This assumption is valid 
in most of the parameter range. Lebedev et al. reported 
spinning motion, where the major axis rotates out of the 
vorticity plane, at very large 7* and large 77*^^ (7* > 2 
and 77i*jj > 8 at = 0.9) [17]. Thus, spinning dynamics 
of microcapsules is also expected at large 77*^^ but it is 
beyond the scope of our present study. 

Here, we considered membrane with constant surface 
area. This condition is valid for lipid vesicles and cells, 
because the lipid membrane is almost incompressible. 
Compared to them, polymer membranes are area com- 
pressible. For microcapsules with a prolate shape, the 
effects of area expansion can be included by a slight mod- 
ification of F and factors for < i < 3 in Eqs. ([J), ([2]), 
and ([5]). For a spherical microcapsule in the absence of 
flow, however, area expansion is essential, since no defor- 
mation from a sphere is allowed without area or volume 
change. The spherical capsules have energy minimum at 
Q!D = so that no saddle point appear even at /i = 0. 
Thus, the spherical capsules likely do not have 'TB with 
TT' phase. 

Capsules with V* < 0.7 may have multiple minima in 
the free-energy landscape; they would show more compli- 
cated phase behavior. Fluid vesicles with V* ~ 0.6 have 
three (meta-)stable states: discocyte, prolate, and stom- 
atocyte [65|. Shear flow induces the elongating shape 
transition from discocyte or stomatocyte to prolate as 
well as the shrinking transition from prolate to discocyte 
[201] . The transition dynamics between discocyte and 
prolate can be well reproduced by two differential equa- 
tions of shape parameter and inclination angle 6 similar 
to Eqs. ([1]) and (j4|). However, the transition dynam- 
ics from stomatocyte to prolate cannot be reproduced 
by the above two equations. At least two shape param- 
eters are necessary to describe the direct shape change 
(no t via discocyte) obtained in our previous simulations 
[20| . When shear elasticity is added to discocyte vesi- 



cles, these metastable states become less stable. In our 
previous paper, we investigated the dynamics of RBCs 
at low shear rate by dynamic equations with one shape 
parameter [4l|. At large shear rate with large 77*^, RBCs 
can become ~ 0, where one shape parameter may not 
be sufficient to describe the RBC shape. 



SUMMARY 



We have systematically studied the dynamics of pro- 
late microcapsules in steady shear flow. The capsules 
exhibit four types of motions, tank-treading, tumbling, 
and two types of swinging oscillations. When the po- 
tential has a saddle point at the phase angle = 7r/2 
(low shear elasticity), these motions can occur simulta- 
neously with or without synchronization. As shear rate 
increases at small viscosity ratio 77;*^, the capsules show 
the transition from tumbling with no tank-treading mem- 
brane rotation to tank-treading mode via coexistence of 
tumbling and tank-treading rotations. Swinging phase 
appears at large 77;*^ and large shear rate 7*, where two 
type of swinging oscillations coexist. When the poten- 
tial has no saddle point (high shear elasticity) , the phase 
diagram is drastically changed. The coexistence phases 
disappear and only two phases (tank-treading and tum- 
bling phases with no tank-treading) remain. 

We have investigated only prolate capsules with a ho- 
mogeneous membrane in this paper. However, our model 
can be applied easily to other capsules by the modifica- 
tion of capsule potential F. Thus, it can cover a wide 
range of soft objects from fluid vesicles to microcapsules 
and cells. 

The coexistence of tumbling and tank-treading rota- 
tions has not yet been observed in experiments and sim- 
ulations. Our study predicts it would occur for low shear- 
elastic capsules or capsules with dimples like RBCs. Fur- 
ther experiments and simulations are expected to confirm 
our predictions. 
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Appendix A: Keller-Skalak Theory 

In the KS theory, a vesicle is assumed to have a fixed 
ellipsoidal shape. 



2 / 0:3^2 
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where denote the semi- axes of the elhpsoid, and the 
coordinate axes Xi point along its principal directions. 
The xi and X2 axes, with ai > 02, are on the vortic- 
ity (xy) plane, and the 2:3 axis is in the vorticity (z) 
direction. The maximum lengths in three directions are 
Li = 2ai, L2 — 2a2, and L3 — 203. The velocity field on 
the membrane is assumed to be 

v-=^(-^X2,^xi,0). (A2) 

Equations ([2]) and ([3]) of vesicle motion are derived from 
the energy balance between supply from the external fluid 
of the vesicle and dissipation inside the vesicle and on the 
membrane. The factors /o, /i, /2, and /s are given by 

/o = 2/(01/02+02/01) = (l-a^)/(l + a^), 

fi = 0.5(01/02 - 02/01) = 2aD/(l - a^), 

/2 - 0.5gial + al), 

h = 0.5£;,i?o/(/iV), 

3= / {al+s)-^/\al+s)-^/\al + s)-^/^ds, 
Jo 

ai = ai/(oia203)^/^. 



where is an integral of shear stress over the membrane 
surface [l^, [3l|. At V* = 0.9, we numerically obtained 
/2 = 0.3185 + O.Gal + 14a^ + exp{124(aD - 0.39)}. 



Appendix B: Estimation of shear elastic energy 

For quasi-spherical capsules, the capsule shape can be 
expanded in spherical harmonics Yi^m as R — Ry{1 + 
J2i m '^i,mYi.m)- Ellipsoidal shapes are described with I = 
0, 2 modes as i? = i?v{l - As/87r + ^2,2(^2,2 + ^2,-2) + 
U2,oY2fi} - i?v{l-As/87r+v/15/87rsin2 {9) COs{2(p)u2,2 + 
(v/5/V/4)(3cos2(6') - 1)m2,o} with u^q + 2u\ ^ = As/2 
in the spherical coordinate representation {R,6,(p). A 
prolate shape with U2.0 — '"2'o ~ —a/ As/2/2 and 1/2,2 — 
^2^2 = -\/3As/ 4 is chosen as a rest state in the absence of 
flow. The displacement is assumed in the radial direction: 

U = UrOr with Ur = {U2,2 - W2^2)(^2,2 + ^2,-2) + ("2,0 - 

W2i))^2,o- The shear elastic energy fi J {si^i — £2,2)^/2 + 
2ei,2^dS is numerically calculated, where Sij is the two- 
dimensional strain tensor on the membrane surface. 
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